Nucleation of bended vortices in Bose-Einstein condensates in elongated traps. 
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We study the generation of vortices in rotating axially elongated magneto-optical traps, a situation 
which has been realized in a recent experiment (K. W. Madison, F. Chevy, W. Wohlleben, J. 
Dalibard, Phys. Rev. Lett. 84 806 (2000)). We predict that at a critical frequency the condensate 
experiences a symmetry breaking and changes from a convex cloud to a state with one bended vortex. 
We also discuss several effects which enlarge the critical frequency with respect to other geometries 
of the trap: these are, (i) the failure of the Thomas-Fermi approximation on the transverse degrees 
of freedom of the condensate, (ii) the enhancement of the transverse asymmetry of the trap by 
means of rotation and (iii) the yet unobserved bending of the vortex lines. 
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Since the creation of the first Bose-Einstein conden- 
sates (BEC) using ultracold dilute gases [Q and the suc- 
cess in describing them using simple but accurate mean 
field theories, it has been a goal of the field to produce 
and understand the Physics behind vortices in these con- 
densates. 

Vortices, or vortex-lines we should rather say, consti- 
tute the most relevant topological defect that one finds 
in Physics. They basically consist on a twist of the phase 
of a wave function around an open line and they are typ- 
ically associated to a rotation of a fluid whatever the 
fluid is made of (real fluids, optical fluids, quantum flu- 
ids, ...) d|. Furthermore, vortices are one of the means by 
which quantum systems may acquire angular momentum 
and thus react to some perturbations of the environment. 
They have already been predicted, observed and studied 
in the superfluid phase of 4 He and are indeed known to be 
the key to some important processes inside these systems, 
such as dissipation, moments of inertia and breakdown 
of superfluidity. This is why extensive research on vortex 
generation, stability and dynamics has been conducted 
in the field of BEC in the last years ||-|| . 

In a recent experiment performed by a group from the 
Ecole Normale Superieure (ENS) (?]] vortices are created 
by rotating an elongated trap slightly deformed in the 
trasverse dimensions. This rotation is performed at a 
constant angular speed, ^, which is the key parameter 
for the appearance of one or more vortices. Many in- 
teresting phenomena have been observed in these exper- 
iments, but we are concerned with some of them which 
are controversial: (i) the first vortex nucleates at a crit- 
ical rotation speed or critical frequency, fii, larger than 
the estimate of Qj F = ^ln^p §, (ii) vortex 

cores seem partially filled and (iii) after the nucleation of 
the first vortex the angular momentum seems to depend 
continuously on the rotation speed. 

In this letter we consider these effects within the frame- 
work of a mean field theory for dilute gaseous condensates 



in rotating traps. In addition to an explanation of the ex- 
perimental features we find the surprising fact that above 
f^i a symmetry breaking changes the ground state (GS) 
into a bended vortex [Fig. |l|. We will discuss this phe- 
nomenon later. 




FIG. 1. Shape of the ground state of the Hamiltonian (g) 
for Q = 0.75cj ± ,7 = 350, Ng = 1000 and e = 0. (a) Trasverse 
section, (b) Density plot. One quarter of the condensate has 
been removed to allow the direct observation of the inner 
structure of the vortex line. 

The model- For most of the current experiments it 
is an accurate approximation to use a zero temperature 
many-body theory of the condensate. In that limit the 
whole condensate is described by a single wavefunction 
^(r, t) ruled out by a Gross-Pitaevskii equation (GPE). 

In the ENS experiment the trap is initially harmonic 
with axial symmetry, but then a laser is applied which de- 
forms it and makes it rotate with uniform angular speed 
Q. By using the mobile reference frame which rotates 
with the trap, we can describe this experiment using a 
modified GPE which reads 
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Here the ENS trapping potential is given by 
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and L z = i [xd y — yd x ) is a linear operator which repre- 
sents the angular momentum along the z-axis. 

In Eq. ([y) we have already applied a convenient 
adimensionalization which uses the harmonic oscillator 
length, a±_ = \/Ti/mRbU±_ and period, r = coj 1 . Using 
these units the nonlinear parameter becomes g = as/a±, 
where as — 5.5 nm is the scattering length for 87 Rb, 
which is the gas used in Ref. @. Following the ex- 
periment we will take uo z = 2ir x 11.6 Hz and co± = 
2tt x 232 Hz. Through the paper we will also use a fixed 
value of Ng = 1000 which corresponds to a few times 10 5 
atoms. For the small transverse deformation of the trap 
we have used e = 0.03, 0.06, 0.09 (e = 0.03 is the closest 
one to the actual experiment). 

The norm, N[ip] =/|^| dr, which can be interpreted 
as the number of bosons in the condensate is a conserved 
quantity as well as the energy 



:A + V (r) 



i/jdr. (3) 



Each stationary solution of the GPE of the form 
Vv( r 5 1) = e _z/i V( r )? is a critical point of the energy sub- 



N 



ject to the constraint of the norm, ^ 

Ground states.- It is physically plausible to assume 
that the ENS experiment actually produces condensates 
which are close to the absolute minimum of the energy 
subject to some constraints. In practice, this means that 
any theoretical effort to reproduce the experimental re- 
sults must be based on a method to find these ground 
states of Eq. (|J) for given norm and rotation speed. To 
achieve this goal and to be able to make quantitative 
predictions we have worked numerically with the energy 
functional [Eq. (||)], using a method known as Sobolev 
gradients. This method is a preconditioned descent tech- 
nique which allows to get accurate approximations to the 
ground state where other minimization procedures fail. 
We have applied this method over a Fourier basis with 
32 x 32 x 64 modes, which is enough to obtain accurate 
estimates of macroscopic values, such as the angular mo- 
mentum and the lowest critical frequencies. Nevertheless, 
the most controversial results have been confirmed on a 
grid with 64 x 64 x 128 modes. 

Our simulations are summarized in Figs. |l| and ^|. The 
first striking prediction is that vortex lines are nucleated 
with a deformed shape, like a tight string which is pulled 
from the middle as rotation is increased [Fig. [IJ. The 
second important result [Fig. ||(b),(d)] is that even for 
a weak trap anisotropy, e = 0.1 [Fig. ^(d)] a transverse 
section of the condensate shows a very deformed cloud: 
the rotating condensate amplifies the small transverse de- 
formation which is used to transfer angular momentum 
to the condensate. 



Furthermore, Fig. || reveals that our numerical simula- 
tions confirm several other features already found in the 
experiment. First, the critical frequency for the nucle- 
ation of a vortex Qi is larger than what it is expected in 
the spherically symmetric trap. According to Fig. ||(a) 
in our simulations (see parameters above) the condensate 
remains in a vortexless state until Q = 0.7(1) x cu_\_/27r, 
that is, fii c± 162 Hz, which is close to the value of 147 
Hz obtained experimentally [Q , but far from the Thomas- 
Fermi estimate for radially symmetric clouds. 
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FIG. 2 

tor, a = (y 2 )/(x 2 ), of the GS of Ety) [Eq. ([J] as a func- 
tion of Q. Solid, dashed and dotted lines correspond to 
£ = 0.03, 0.06, 0.09 respectively, (c-d) Transverse density 
plots of a GS without vortex (Q = 0.65, £ = 0.1) and a GS 
with one vortex (Q = 0.75, £ = 0.1). Pictures are 12x12 
adimensional units large. 

The second remarkable feature is that, at Q = the 
angular momentum per particle (L z ) /N jumps to a value 
which is smaller than 1 and which grows continuously 
until a second vortex is nucleated [Fig. ||(a)]. 

Thus, the numerical simulations are in good agreement 
with the experience and both contradict previous ideas: 
(i) the nucleation frequency is much larger than expected 
in Thomas-Fermi theories, and (ii) the angular momen- 
tum has a continuous, nonlinear dependence on Q. To 
obtain more insight on the physical reasons for this be- 
haviors we will analyze several simplified situations. 

Critical frequencies for axially symmetric traps .- Let 
us first consider the case of an axially symmetric trap. 
For these traps it is usually thought that at a certain Vt = 
Qi the system experiences a transition from a vortexless 
state to a symmetric state with one vortex. Both states 
do also exist in the nonrotating case and have energies 
and E\ > Eo, respectively. This leads to the usual 
estimate of the critical frequency Q\ c± E\ — Eq. 
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FIG. 3. (a) Critical frequency, Q± for a spherical trap 
(7 — — 1? solid line), a cigar shape trap (7 = 19, upper 

dashed line) and a pancake trap (7 = 1/19, lower dashed 
line), (b) Critical frequency Qi (solid line) and eigenvalues 
arising from the linearization of the energy (circles) for a trap 
with 7 = 19. It is remarkable that there is always a negative 
eigenvalue which is not suppressed for Q = Qi. (c) Vortex, 
(d) main destabilizing mode and (e) bended vortex line arising 
from a combination of both, for a trap with 7 = 19. 

Taking the previous idea for granted and to learn how 
the elongation of the trap affects the nucleation of the 
first vortex, we have searched the vortexless state and 
the state with one vortex in different axially symmet- 
ric traps (e = 0), computing {Eq : Ei, Oi}. The result 
of these calculations is shown in Fig. ||(a) for different 
elongations. There we see that the critical angular speed 
of a cigar shape trap is larger than that of spherically 
symmetric and pancake traps. The explanation is sim- 
ple: in an elongated trap, uo z <^i, the system expands 
preferably along the axial dimension in order to accom- 
modate more bosons and counteract the interaction en- 
ergy. As uo z is made smaller this tendency becomes more 
evident, since the transverse shape no longer satisfies the 
Thomas-Fermi approximation in which the second order 
derivatives are negligible. Instead, the transverse shape 
of the cloud is closer to a harmonic oscillator wavefunc- 
tion, and as such the energy required to introduce a single 
vortex, fii, becomes larger than in the spherically sym- 
metric case. In other words, the more elongated the trap 
is, the more dilute the condensate seems on the XY plane, 
the closer Qi becomes to its linear limit uu±, and the faster 
the condensate must spin to accept a vortex. This effect, 
which can be characterized as a dynamical separation of 
longitudinal and transverse degrees of freedom || , mani- 
fests in a departure of Qi from any Thomas-Fermi based 
results, which correspond to a regime ruled by the trans- 
verse self-interaction of the condensate. 

Bended vortices in axially symmetric traps.- An impor- 
tant feature of elongated traps is that they require very 



little energy to produce longitudinal excitations of the 
condensate ||. This makes feasible for small perturba- 
tions to induce longitudinal modes and distortions of the 
vortex lines pJ]||L2| . Thus our next step is to perform a 
linear stability analysis of the straight vortex line which 
we have already found. 

We are looking for a quadratic expansion of the energy 
of a perturbed vortex. Since the vortex Vw* is a sta- 
tionary state, this energy must read, up to second order 
perturbations 

E(t/j vor t + e5)~E 1 -n + e 2 Y^ (Km - mil) \c ijfn \ 2 (4) 

i 

Here, A; ?m are the different eigenvalues which arise from 
linearizing the GPE [Eq. ([l])] around tp vor t, rn are their 
vorticities with respect to t/j vor t (See ||), and c^ m is the 
weight of the \i,m) mode in the expansion of S(r). By 
finding the lowest negative eigenvalue, A^_i < 0, we 
know that at a rotation speed ft m eta = the en- 

ergy (|J) becomes positive for all perturbations and the 
vortex is locally stable or metastable. 

These calculations are summarized in Fig. ||(b). There 
we plot all eigenvalues — A^_i and find that there is a 
number of them such that — A > {E\ — Eq). This result 
suggests that there exists a range of speeds Q G (fii = 
Ei — Eo,Q m ) in which the system nucleates a vortex, 
but that vortex does not have a symmetric shape. In- 
stead, the stability analysis points out that this vortex 
is bended by its ends, forming a curved vortex line [Fig. 
||(e)]. Although this is a perturbative prediction arising 
from the linear stability analysis, these solutions are con- 
sistent to what we formerly obtained by minimizing the 
full Hamiltonian and which was explained above [Fig. 

The existence of a ground state with such a bended 
shape is one of the main predictions of the paper, since it 
presents a symmetry breaking due solely to rotation. In 
some sense the symmetry is broken twice: First because 
a phase asymmetric structure, the vortex, appears, and 
a second time because the vortex bends in some specific 
direction. We have also observed that as the frequency 
is increased from Qi the vortex bending becomes less 
and less evident (And (L z ) approaches TV). However, 
when the threshold ^2 for the nucleation of two vortices 
is reached the bending becomes again very important. 

The existence of bending in the vortex lines may be 
checked in the experiment and is probably related to the 
experimental observation that vortex lines appear to be 
partially filled when viewed from above. 

Role of transverse asymmetries.- In Ref. the trap is 
distorted over the transverse directions to induce some 
mechanical response of the condensate, an effect which 
must be taken into account. 

As it was already shown [|ll| any initial asymmetry of a 
condensate is strongly emphasized by the rotation of the 
trap, an effect that may be estimated analytically in the 
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linear and Thomas- Fermi limits. These estimates lead to 
a value of the angular momentum for a vortexless conden- 
sate and a different one for a condensate with a single, 
centered, unit-charge vortex line. The main difference 
with respect to the symmetric case is that a vortexless 
state may acquire some angular momentum and that the 
nucleation of a vortex line gives an increase of the total 
angular momentum which is inversely proportional to the 
asymmetry factor, e. From a practical point of view this 
means that the system may "decide" to wait for a larger 
rotation speed > fii jS2/m before nucleating a vortex. 
Actually what happens is that fii/min(a; x , u y ) increases 
with respect to e, but decreases [Fig. g|. 

Although it is clear that the effect of transverse asym- 
metries is to provide a larger critical frequency than ex- 
pected it is not possible to use any analytic estimate of 
the critical speeds because of the elongation of the trap 
and its transverse deformation. Furthermore, due to the 
bending of the vortex lines it is not possible to work with 
the corrections from |ll[ to get a reliable estimate of Qi. 

Discussion and conclusions.- We have used the GPE to 
study theoretically the nucleation of vortices in an elon- 
gated trap similar to that of the ENS experiment |7j . We 
have found several mechanisms which lead to an increase 
of the angular speed which is required to create those vor- 
tices and have found the remarkable feature that vortex 
lines are bended by their ends. This provides some kind 
of symmetry breaking bifurcation in the case of a sin- 
gle vortex |l3| . Recalling previous experience with liquid 
helium, we suspect that this bending has some impor- 
tant consequences, and that it might be responsible for 
the breakdown of superfluidity in BEC by a mechanism 
which consists in several vortex lines joining together and 
forming a turbulent pattern. 

Using only the special geometry of the experiment we 
have obtained a value of fti = 0.7(0)cjx which is close but 
above the experimental value. However, this means that 
the experiment implies a systematic negative shift both 
in Qi and in the angular speed at which the superfluid 
escapes, Q c . In the experiment tt c ~ 2tt x 210 Hz, while 
the theory predicts Q c = uo_\_(l — e). If we trust the value 
of Q c from the experiments Q c and make cu_\_ = 2tt x 210 
Hz, our prediction is Q\ = 2tt x 147 Hz, which coincides 
exactly which the experimental measures j?J . 

The fact that simple mean field theories are still valid 
when vortices enter into the dynamics, and that we may 
still expect quantitative predictions based on this simple 
models is interesting and has been shown in many other 
phenomena The striking result that the combined 
effect of asymmetry and nonlinearity may lead to the 
phenomenology here studied is another example of the 
richness of behaviors that one may expect from these 
nonlinear Quantum mesoscopic systems. 
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